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Equivalent representations of higher-order Hamiltonian guiding-center theory 
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Two complementary representations of higher-order guiding-center theory are presented, which 
are distinguished by whether higher-order corrections due to magnetic-field nonuniformity appear 
in the guiding-center Poisson bracket or the guiding-center Hamiltonian. The equivalence of these 
two representations implies that the guiding-center gyroaction adiabatic invariant and the guiding- 
center Hamilton equations of motion are representation-invariants. By using a new perturbative 
hierarchy in guiding-center theory, second-order corrections are introduced in Hamiltonian guiding- 
center theory without the need to calculate the guiding-center transformation at second order. 



Guiding-center theory [H, H[ is a powerful paradigm for 
the reduced description of charged-particle dynamics in 
weakly-nonuniform magnetic fields, which yields a deep 
understanding of the physics of magnetic confinement in 
laboratory and space plasmas over long time scales 
In his pioneering work, Littlejohn [ij used Lie-transform 
perturbation methods § to derive noncanonical Hamil- 
tonian guiding-center equations of motion that contain 
first-order corrections associated with the magnetic-field 
nonuniformity. These first-order corrections not only 
have theoretical importance in constructing a hierarchy 
of adiabatic invariants (e.g., (6() but also play a signifi- 
cant role in ensuring the accuracy of guiding-center nu- 
merical codes 0] in following magnetically-trapped par- 
ticles in complex magnetic fields. 

The first attempts at keeping second-order corrections 
in the expressions for the adiabatic invariants appeared in 
non-Hamiltonian guiding-center equations 0, |9[ . Hamil- 
tonian guiding-center equations that include second- 
order corrections (appearing in the guiding-center Hamil- 
tonian) were recently derived by Parra and Calvo [ic| |. 
who pointed out that these higher-order corrections were 
required in order to construct a consistent Hamiltonian 
gyrocenter theory [Tl1 | that retains all second-order cor- 
rections due to magnetic nonuniformity as well as fluctu- 
ation amplitude and space-time scales. Because of their 
potential relevance in numerical gyrokinetic applications 
[1 21 ] , a complementary higher-order Hamiltonian guiding- 
center theory is constructed here, which includes second- 
order corrections (appearing in the guiding-center Pois- 
son bracket), by using a new perturbative hierarchy in 
guiding-center theory. 

Modern guiding-center Hamiltonian theory [2j is de- 
scribed in terms of noncanonical guiding-center phase- 
space coordinates z a = (X.,p\\ , H, 9), which include the 
guiding-center position X and the guiding-center paral- 
lel momentum , while the gyro-action J = (mc/e) [i 
(defined in terms of the magnetic moment \i for a par- 
ticle of mass m and charge e) is canonically-conjugate 
to the ignorable gyroangle 9. The guiding-center phase- 



space coordinates z a are constructed by a near-identity 
transformation from local particle phase-space coordi- 
nates Zq = (x,p|| , /xo, 6*o), generated by the phase-space 
vector fields G„ (n = 1,2,...) @: 



z^ + eG^+e' G% 



\ Gi • dG? 



. (1) 



where the ordering parameter e appears through the 
charge renormalization [l[ e — > e~ 1 e. 

Hamiltonian guiding-center theory 0, 0] is based on 
the guiding-center phase-space Lagrangian Q 



T- c x Ao 



dS 



- gc 



H„ r dt, 



(2) 



where the guiding-center push-forward operator T~ c is 
generated by (Gi,G2,---) and the gauge-function S = 
e Si + e 2 ^ + • • • is chosen to simplify the form of the 
guiding-center symplectic one-form r gc . The local phase- 
space Lagrangian Ao = Fo — Hq dt in Eq. {2} is ex- 
pressed in terms of the local particle symplectic one-form 
r = (e _1 e A/c + piiQ b + pxo) - dx, with the magnetic 

field defined asB = VxA = £?b, and the local particle 
Hamiltonian Hq = p| /2m + /iq B. 

The results of the Lie-transform perturbation analysis 
@ are summarized as follows [13|. First, the generic 
form of the guiding-center symplectic one-form is 

F gc ee ^A + n||b) -dX + J (d6 R-dX 



EE - A* 



• dX + J 66, 



(3) 



where the vector field R(X) ensures that the one-form 
60 — R-dX is gyrogauge- invariant In Eq. the 
symplectic parallel momentum 



H, 



P\\ 



(1) 



(4) 



is defined in terms of the higher-order corrections 



n 



(i) 



= — (D 
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b,(6) 
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where (• • • ) denotes gyroangle averaging, p = — G* = 
(cb/eB) X pj^ denotes the lowest-order gyroradius, the 
first-order operator D\ (C) = Gi • dC + VC • p Q is defined 
for an arbitrary vector field C, and we used D n {p\\ b) • b = 

Gn l — p\\ G* • k for n > 1, where n = b • Vb denotes the 
magnetic curvature. 

Next, the guiding-center Hamiltonian is defined as 



H gc = T gc H - — 



where the guiding-center ponderomotive potential 



e 2 * 2 



is defined in terms of the higher-order corrections 



*2 



m 

\ (G!)-d*i, 



(7) 
(8) 
(9) 
(10) 



where we used (p ) = in Eq. (j9]). In Eqs. (J6j) and 
(|g])-(|T0j). the components (G*,G£,G^) are determined 
from the perpendicular spatial and momentum compo- 
nents of the guiding-center symplectic one-form. 

In the perturbation analysis leading to Eqs. ([3]) and 
(0, the generating vector field G„ is chosen so that the 
corrections nj,"' and ^> n are both independent of the 
guiding-center gyroangle 6, which yields a guiding-center 
magnetic- moment adiabatic invariant up to e n+1 . Once 
the guiding-center phase-space Lagrangian @ has been 
constructed to any desired order (we now replace e/e with 
e), the guiding-center Hamilton equations of motion are 
expressed in terms of the guiding-center Hamiltonian ([7]) 
and the guiding-center Poisson bracket constructed from 
the guiding-center symplectic one-form ©: 



{F, G} g 



dFdG _ d£dG 

~d6~dJ ~ ~dJ ~de 

{dpj B 
cb 



BG OF 



dp« dp 



V*G 



eB 



(V*F x V*G) , 



(11) 



where the symplectic magnetic field B* = V X A* is 

c 



B* 



B 



- v x (n,, b 



JR 



(12) 



with BTi = B* • b, while the gyrogauge-invariant gradient 

V* = V + R* d/d6 is defined with R* = R b 8TL\\/dJ. 

The guiding-center equations of motion are now ex- 
pressed in Hamiltonian form as d gc z a /dt = {z a , H gc } gc , 
where the guiding-center operator d gc /dt is defined in 
terms of the particle operator d/dt and the guiding-center 



push-forward and pull-back operators T^ 1 and T g 
erated by (Gi, G 2 , • • • ) as 



dgc T _i ( d 
dt ~ gc \dt g 



The reduced guiding-center equations of motion 
'dH\\ 

II / " ■ -| 



gen- 



(13) 



d gc ~X- 



dt 



dp 



1 dH gc B* cb ^ TT 

^pJb; + ^ xVH ^ (14) 



'ii 



/an,, v 1 b* 



dt 



V dp 



B 



(15) 



are decoupled from the gyro-motion equations d gc 9 / dt — 



fl + R* • d gc X./dt and d gc J/dt 



dH gc /d9 = 0. 



Moreover, in contrast to non-Hamiltonian guiding-center 
equations (e.g., [l[) the reduced guiding-center Hamilton 
equations (fl4l) - (IT5l conserve energy (d gc H gc /dt = 0) and 
satisfy the guiding-center Liouville theorem 



d gc X 
dt 



d 

dp\\ 



d g cP\\ 
dt 



= 0, (16) 



where J gc = BT, (dH\\/dp\\) is the Jacobian for the 
guiding-center phase-space transformation (p}. Note 
that when II || and ^ include nth-order correction terms, 
the guiding-center Hamilton equations (|Ti)) - (fT5"|) contain 
terms of order n + 1 through B* and V^. 

While the reduced guiding-center Hamilton equations 
(TT^1l - ([T5|l are completely general and are valid at all or- 
ders in e, they can be simplified if the higher-order cor- 
rections are placed either entirely in the guiding-center 
Poisson bracket or in the guiding-center Hamiltonian. 
In the former case, called the symplectic representation, 
the symplectic parallel momentum (|4]) contains all the 
higher-order corrections (i.e., Tl^ for n > 1), while 
the ponderomotive potential © is simply ^> = fj,B (i.e., 
"fn = for n > 1). In the latter case, called the Hamilto- 
nian representation, the symplectic parallel momentum 
IT || replaces p|| as a dynamical variable (i.e., nj"' 1 = for 
n > 1), while the ponderomotive potential © contains 
the higher-order corrections (i.e., ^ n ^ for n > 1). 

These two complementary representations are gener- 
ated by two different phase-space transformations ([1]) 
leading to the guiding-center coordinates (X,pii,|U, 6) 
generated by (Gi, G2, ■ ■ • ) in the symplectic representa- 
tion or to the guiding-center coordinates (X, py , ~p, 9) gen- 
erated by (Gi, G2, • • • ) in the Hamiltonian representation. 
These representations are said to be equivalent if 



( X ,P||,M 



(x,nN,u, 



(17) 



so that the Jacobian for the guiding-center transforma- 
tion leading to the Hamiltonian representation is J gc = 
Bn, since dHn/dpu = 1. The representation equivalence 



3 



(fT7| implies that, for a ^ p\\, we have the following rela- 
tions between generating vector fields 



G a 2 



GO. 
1 ) 



1(D OGl 
2 U ll dp,, 



while, using Eq. (|4]), we have 



rfp\\ — r p « 

\_r^ = 1 

7=7Pii — r iP 



n 



(1) 



(18) 
(19) 

(20) 



n 



(2) 



The representation-invariance of the magnetic moment 
Ji = \i implies that, once the corrections Ilj, are found 
in the symplectic representation, the corrections $ n are 
immediately found in the Hamiltonian representation as 



* 2 



m 

P\\ 



(i) 



n 



(2) 



1 

2m 



n 



(i: 



(22) 
(23) 



Hence, once a higher-order Hamiltonian guiding-center 
theory is derived in one representation, one can also 
construct a theory in the complementary representation. 
More importantly, it can also be shown 13[ that the 
guiding-center Hamilton equations (jT4j)- (JT5J are identi- 
cal in all equivalent representations. 

We now proceed with the construction of a higher- 
order Hamiltonian guiding-center theory that retains up 
to second-order corrections. While it appears that the 

(n) 

derivation of the nth-order corrections II and re- 
quires obtaining the guiding-center phase-space transfor- 
mation (fTJ) up to nth order, we now show that, by intro- 
ducing the perturbative guiding-center Banos hierarchy 
in the symplectic representation, we obtain the remark- 
able fact that the nth-order Hamiltonian structure can be 
determined from the guiding-center phase-space transfor- 
mation (fTJ) up to (n — l)th order only. 

The guiding-center Banos hierarchy is constructed in 
the symplectic representation as follows. Using the func- 
tional definition (|13p for d gc /dt, we introduce the identity 



dx 
~dt 



(J-1 \ 



dt 



dt 



expressed in terms of the guiding-center velocity d gc X./dt 
and the guiding-center displacement velocity d gc p gc /dt 
(which includes the polarization velocity d(p gc ) /dt). 
Here, the guiding-center displacement is expanded as 

Pgc = Tgc 1 * x = epo + e2 Pi + g3 P 2 + • • • . ( 2 5) 
where the higher-order corrections 

Pi = - GJ - \ Gi ■ d Po , (26) 
p 2 = - Gf - G 2 ■ d Po + i d ■ d(Gi • d Po ) (27) 



satisfy (p n ) =/= and p n • b ^ 0. By taking the dot prod- 
uct of Eq. (|24|) with m b (defined in guiding-center phase 
space), we obtain the guiding-center Banos hierarchy 



T„ 



b • m — - 

dt 



m b • 



d gc X 
dt 



d gc p g 
dt 



(28) 

where we used the identity / (T gc \g) = T g ^(T gc f g). 
The left side of the guiding-center Banos hierarchy (|28|) 
is expanded in powers of e as 



t: 



eG^ 



V\\ 



-ep b .Vb-p J . + -.- ,(29) 
where the gyroangle-averaged first-order correction 

; Po .vb. P± \ = j T (30) 



represents the Banos parallel momentum [{J [TJ], de- 
fined in terms of the twist of the magnetic-field lines 
t = b • V X b. The right side of the guiding-center Banos 
hierarchy 11280 involves 



nib • 



dt 



dt 



= P\\ - eP|| 
= e mb • I fl 



dU \\ ] 



dp\\ 

dpi d p 
86 dt 



(31) 
,(32) 



where b • dp /d9 = and the guiding-center operator 
(fT3"|) is expanded as 



dt 



i:)0 



dt 



dt 



(33) 



The guiding-center Banos hierarchy (|28p is a perturbative 
hierarchy that relates (Gn ) with c?n,| /dp\\ at each order 
n > 1. At first order, the guiding-center Banos hierarchy 
yields the gyroangle-averaged relation 



<G?> = 



P\\ 



an i| 1} 
dp\\ 



Jr. 



(24) When it is combined with Eq. ©, with Eq. ((30)) : 



we obtain the first-order Banos equation 

4 ("" n " ') " - 5 



Jr, 



(34) 



(35) 



(36) 



where the unknown (G^ 11 } has been eliminated. Inte- 
grating Eq. (|36l) with respect to pn yields the first-order 
correction 



n 



(i) 



Jr. 



(37) 
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This solution means that, since (G^ 11 } = J r in the sym- 
plectic representation, the first-order Bahos parallel mo- 
mentum J r is hidden in the definition of the guiding- 
center parallel momentum. Note that the first-order sym- 
plectic momentum (I37|) was obtained without knowing 
the details of the first-order perturbation analysis [H, 0] • 
At second order, we similarly obtain the second-order 
Banos equation 



d 
d P]] 



(piinj 2) ) = 2pnej||K| 2 - ft, 



(38) 



where g\\ = p^/(mfl) and the scalar field 



3 k ■ 



V- 



- VlnS 
K + b (V-b) 



l r 

4 



(V-b) J 



(39) 



where qh = 11^/ (mil) and R* = R + \ (r + Q\\ P2) b, 
while Eq. (fT4"|) involves the factor 



dlL, 

Spy 



= 1 



(43) 



where the first-order corrections vanish. 

In the Hamiltonian representation (where II|| acts as a 
dynamical variable), on the other hand, the equivalence 
relations (j2"2")) - (l2T>l) are used to construct the guiding- 
center pondcromotive potential 



n? 



1 



1 Jt 5 



8 rafi 



2m 



I (JP2 
in 



4 n 2 



(44) 



involves terms of second order in magnetic-field nonuni- 
formity. Integrating Eq. (|38[) with respect to yields 
the second-order correction 



n 



(2) 



p\\ q\ \ k ? 



J Q\\ ft) 



(40) 



which is obtained solely from the details of the first-order 
perturbation analysis [2], |4| . 

In the symplectic representation (with "J = fiB), the 
guiding-center symplectic parallel momentum is thus 



n,. 



p\\ 



i4w 



Q\\ & 



(41) 



so that, up to second order, the guiding-center symplectic 
one-form is 



6 -A 

c 



1 



Q\\ \ K \ 



•dX 



J (66 - R* • dX) , 



(42) 



All second-order terms in Eq. (|44|) appear in Eq. (135) of 
Parra and Calvo [10(. However, since they use a mixed 
representation (i.e., TVtP ^ and ^ 2 7^ 0), which may 
not be equivalent to ours, there are additional terms as 
well (which do not satisfy the Banos hierarchy). 

In summary, we have presented two equivalent rep- 
resentations of higher-order Hamiltonian guiding-center 
theory that retains second-order corrections and satisfies 
the perturbative Banos hierarchy ([28]) derived from the 
guiding-center identity (I24|) . These representations yield 
identical Hamilton guiding-center equations, which can 
be used to construct a consistent Hamiltonian gyrocenter 
theory used to study the long-time magnetic confinement 
of charged particles in the presence of electromagnetic 
fluctuations. 
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